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A Carbon Isotope Mass Balance for an Anoxic Marine Sediment: Isotopic Signatures 

of Diagenesis 

ABSTRACT 

A carbon isotope mass balance was determined for the sediments of Cape 
Lookout Bight, NC to constrain the carbon budgets published previously (Martens and 
Klump, 1984; Martens et al., 1993). The diffusive, ebullitive and burial fluxes of 
SC0 2 and CH 4 , as well as the carbon isotope signatures of these fluxes, were 
measured. The flux- weighted isotopic signature of the remineralized carbon (-18.9 ± 
2.7 per mil) agreed with the isotopic composition of the remineralized organic carbon 
determined from the particulate organic carbon (POC) fi* 3 C profiles (-19.2 ± 0.2), 
verifying the flux and isotopic signature estimates. 

The measured fi 13 C values of the SC0 2 and CH 4 diffusive fluxes were 
significantly different from those calculated from porewater gradients. The differences 
appear to be influenced by methane oxidation at the sediment-water interface, although 
other potential processes cannot be excluded. 

The isotope mass balance provides important information concerning the 
locations of potential diagenetic isotope effects. Specifically, the absence of downcore 
change in the 6 13 C value of the POC fraction and the identical isotopic composition of 
the POC and the products of remineralization indicate that no isotopic fractionation is 
expressed during the initial breakdown of the POC, despite its isotopically 
heterogeneous composition. 
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INTRODUCTION 


Continental shelf and slope sediments are the dominant sites of carbon cycling 
in marine sediments, comprising only 10% of the seafloor but accounting for 
approximately 90% of the organic carbon remineralization (Henrichs and Reeburgh, 
1987). These environments are thus critical to the global carbon cycle. Carbon 
remineralization rates in continental margin environments, however, are often difficult 
to quantify due to the temporal and spatial variability in processes. 

Organic carbon reaching the sediments can be remineralized via a succession of 
oxidative processes requiring oxygen, nitrate, iron, manganese and sulfate as electron 
acceptors (Mechalas, 1974). After these oxidants are depleted, if sufficient labile 
organic carbon remains, it may be fermented to CH 4 and CO 2 (Claypool and Kaplan, 
1974). A common approach to estimating the rate of organic matter decomposition is 
the measurement of the sediment water fluxes of the remineralized products, SCO 2 
and CH 4 and/or the oxidants (Henrichs and Farrington, 1984; Berelson et al., 1987; 
Mackin and Swider, 1989; McNichol et al., 1991; Reimers et al., 1992; and others). 
Significant uncertainties in estimates of sediment-water carbon fluxes are 
by temporal variability in the input of organic matter and associated 
rcmir p ra b 7 atin n rates at individual sites (Martens and Klump, 1984; McNichol et al, 
1991; Reimers et al., 1992). One approach to address these uncertainties is to 
determine a carbon mass balance of the various sources and sinks at individual field 
sites on time scales appropriate to the site (Martens and Klump, 1984; Berelson et al, 
1987; Martens et al., 1993). In temporally variable systems, however, carbon budgets 


2 


may not be easily constrained without numerous measurements of fluxes. 

An isotope mass balance, accomplished by measurement of the 6 C of the 
individual carbon fluxes and determining the flux-weighted averages of each of the 
carbon pools, can serve to constrain further the fluxes and rates. Carbon isotope mass 
balances have been attempted in lake productivity studies (Quay et al., 1986; Herczeg, 
1988) and in estuarine studies to quantify seasonal organic matter fluxes to the estuary 
floor (Lucotte et al, 1991). Alperin (1988) determined an isotope mass balance for the 
sediments of Skan Bay, Alaska in order to identify the relative rates of degradation of 
different organic sources. On a global scale, an isotope mass balance of atmospheric 
methane has been used to estimate source inputs (e.g. Craig et al., 1988; Quay et al., 
1991). Most recently, Quay et al. (1992) and Tans et al. (1993) have attempted to use 
SC0 2 6 13 C measurements of surficial seawater to estimate increases in oceanic 
uptake. All of these isotope studies have attempted to use a mass balance to estimate a 
flux or rate that was not easily measured. 

Our approach has been to use isotope measurements of all of the major carbon 
reservoirs and fluxes of the sediments of Cape Lookout Bight, NC, to verify the 
carbon budget proposed for this site (Martens and Klump, 1984; Martens et al., 

1993). Measurement of the isotopic composition of the individual fluxes and 
reservoirs also allows us to address specific questions about processes at the sediment- 
water interface such as methane oxidation and to characterize diagenetic isotope effects 
caused by selective degradation. Cape Lookout was chosen because carbon mass 
balances for this site have been published previously (Martens and Klump, 1984; 
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Martens et al., 1993) and remineralization rates vary seasonally in a predictable 
fashion (Martens et al., 1986). 

Remineralization processes (e.g. microbial methane production) can alter the 
isotopic composition of the individual carbon fluxes, causing a measurable isotope 
fractionation between the CH 4 and the C0 2 . Diagenetic isotope effects associated with 
selective remineralization of isotopically distinct organic fractions have also been 
hypothesized to alter the isotopic signature of the fluxes of the remineralized fraction, 
and by mass balance, the buried carbon fraction as well (Spiker and Hatcher, 1984; 
1987; Alperin, 1988; Fischer, 1989; Benner et al., 1991). In principle, the 12 C and 
13 C mass balance of all identified carbon fluxes (CH 4 , SC0 2 , corrected for carbonate 
dissolution and precipitation) should be the same as the isotopic composition of the 
remineralized carbon fraction determined from the particulate organic carbon isotope 
profiles. 

FIELD SITE 

Cape Lookout Bight (CLB) is a small (1 km 2 ) back barrier island lagoon with a 
water depth of 7 meters at the deepest location (Martens and Klump, 1980). All 
samples were collected at Station A-l, where previous carbon mass balance has been 
established (Martens and Klump, 1984; Martens et al., 1993). Circulation within the 
Bight is controlled predominantly by tidal flow (Martens and Klump, 1980; Wells, 
1988). The water column above station A-l remains oxygenated all year (Bartlett, 
1981). 

Cape Lookout Bight acts as a trap for fine-grained sediments and organic 
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debris transported seasonally from the Atlantic shelf during storms (Canuel et al., 

1990) and from the shallow back barrier lagoons of coastal North Carolina by ebb 
tidal flows (Martens and Klump, 1980; Chanton et al., 1983). Sediment in Station 
A-l accumulates at a rate of 8.4 to 11.8 cm/yr averaged over the upper meter of 
sediment (Chanton et al., 1983; Canuel et al., 1990). Organic carbon content in the 
upper meter ranges from approximately 3 to 4 wt. % (Martens and Klump, 1984; 
Haddad and Martens, 1987; this study) and appears to be derived primarily from 
phytoplankton and seagrass debris (Haddad and Martens, 1987). Bioirrigation and 
sediment mixing by animals are not important processes for most of the year at Station 
A-l (Bartlett, 1982). 

Sulfate reduction is the dominant remineralization process within the upper 10 to 
15 cm of the sediment (Martens and Klump, 1984; Crill and Martens, 1987). 

Methane production occurs below this zone. During summer months (May through 
October) methane ebullition occurs during low tide conditions when hydrostatic 
pressure is decreased (Martens and Klump, 1980). Sediment temperatures can vary by 
more than 20°C through the year and are the driving force behind seasonal changes in 
rates of microbial processes such as sulfate reduction and methanogenesis (Martens 
and Klump, 1984). 

Based on grain size distributions and annually reproducible pore water nutrient 
profiles (Klump and Martens, 1981), Martens and Klump, (1984) suggested that 
"quasi steady state" conditions have existed at Station A-l at least since the early 
1970’s. Haddad and Martens (1987) used lignin oxidation product analysis and 



isotopic measurements of sediment organic matter as well as degradation rates of 
metabolizable materials to document yearly steady state input of organic matter on an 
annual time scale at Station A-l. The agreement in the sedimentation rate determined 
using 210 Pb and 7 Be radiotracers in studies conducted almost 10 years apart (Chanton 
et al., 1983; Canuel et al., 1990) is further evidence of annual steady state sediment 
accumulation at Cape Lookout Bight. Laminations seen in sediment cores in the upper 
5 cm of winter sediment cores and 7 Be profiles (Canuel et al., 1990) indicate, 
however, that sedimentation is not uniform on time scales of less than one year. Pore 
water depth-concentration profiles, fluxes across the sediment-water interface and 
isotope profiles change rapidly in response to increasing temperature, suggesting that 
steady state conditions may not be attained on time scales of weeks. Nevertheless, 
pore water depth integrated and measured rates suggests that the system is near steady 
state on time scales less than a year. 

METHODS 

Cores were collected by divers from Station A-l at 4-8 week intervals from 
February 1986 through February 1987. Samples also were obtained with a Soutar box 
core in March, 1986 and 1988. Porewater samples were isolated with a sediment 

membrane filtration squeezer (Reeburgh, 1967). 

SC0 2 analyses were performed by injecting one to two ml of porewater into 
120 ml evacuated serum bottles capped with crimped, lightly greased rubber stoppers 
(Alltech #6633). Duplicate or triplicate cores were analyzed. In May, 1986, 
duplicate samples from the same core as well as samples from a second core were 
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analyzed. Pore water samples were frozen immediately after collection and 
maintained frozen until analysis (-18 to -56°C). One ml of 1M phosphoric acid 
saturated with copper sulfate was added to each porewater sample immediately before 
analysis. The copper sulfate was added to the sample to precipitate sulfides in order 
to avoid interference from H 2 S. The CO 2 was removed from the bottle by vacuum 
distillation, purified, and collected cryogenically for isotopic analysis. Sample size 
was determined manometrically. Bicarbonate and tank CO 2 standards were processed 
to verify that the procedure was not causing isotopic fractionation and to insure 
complete gas collection. The precision and accuracy of concentration measurements 
were 0.4mM and O.lmM, respectively, based on bicarbonate standards. 

Reproducibility of the triplicate samples from May, 1986 was 0.5mM. 

After purification the CO 2 was collected and sealed in borosilicate tubing and 
its carbon isotopic composition was measured in one of three laboratories: the Stable 
Isotope Laboratory, North Carolina State University; NASA-Ames, California; and the 
Center for Applied Isotope Studies, Athens, Georgia. Interlab comparison among 
these three labs gave comparable results (±0.2 per mil, Blair and Carter, 1992). The 
accuracy and precision of the ECC^ measurements were 0.4 and 0.2 per mil, 
respectively. Reproducibility from triplicate standard analyses was 0.95 per mil. 

Porewater samples were collected on four field dates for analysis of calcium. 
Following acidification and filtration (0.45 iim), samples for calcium were diluted to 
1:20. Samples were analyzed on a Perkin-Elmer 560 Atomic Absorption spectrometer 
(NCSU - Forestry Department Analytical Laboratory) using a 0.5% Lanthanum 
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solution. Duplicate analyses had a precision of 0.2mM. A comparison of samples 
analyzed for Ca ++ at SUNY-Stony Brook (P. Rude, Marine Sciences) and NCSU 
agreed to within 0.2mM. 

Sediment cores were collected for CH 4 isotopic analysis. Cores were 
sectioned in 2-3 cm depth intervals and placed in Mason jars containing 10 to 20 ml of 
1-3M NaOH. Headspace gas was removed from the Mason jars by syringe through a 
rubber stopper (Bellco Biotechnology) that had been fitted into the lid. After 
removing water and C0 2 cryogenically, the methane was converted to C0 2 by passing 
it through an 790°C furnace packed with CuO using helium as a carrier gas 
(Matthews and Hayes, 1978). The resultant C0 2 was purified and collected 
cryogenically. The precision and accuracy of the analyses, based on a methane 
standard (Scott Specialty Gas), were both 0.3 per mil. Duplicate core samples 

expressed a maximum difference of 1.6 per mil. 

Bubble samples were collected as described in Martens et al. (1986). One ml 
samples of the bubble gas were removed via syringe from a serum bottle and 
processed in the same way as the sedimentary methane samples. Tank CH 4 standards, 
duplicate samples from individual bottles and duplicate and triplicate bottles were 
analyzed. Precision and accuracy of this procedure are the same as the values given 
for the sedimentary methane procedure. Reproducibility of samples taken from the 
same bottle were 0.02 per mil and reproducibility of analyses from individual bottles 

collected on the same date were 0.2 per mil. 

Concentration and isotopic composition measurements of the particulate organic 
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carbon (POC) were made downcore on sediments collected in August, 1986. Bomb 
combustion techniques are described in Blair et al. (1987). Precision and accuracy of 
this procedure, based on NBS-20 standards, are 0.3 and 0.5 per mil, respectively. 

Particulate inorganic carbon (PIC) was analyzed on a sediment core collected in 
May, 1987. The sediment was rinsed with deionized water three times and freeze 
dried. Ten to 60 mg of sediment was placed in a 120ml serum bottle and sealed with 
a rubber stopper. The bottle was evacuated and one ml of 0.5M phosphoric acid was 
injected and allowed to react overnight. Vacuum line collection procedures were the 
same as for EC0 2 porewater samples. Duplicate analyses of the 0-1 cm depth interval 
gave PIC concentrations within 0.1% and fi^C values within 0.01 per mil. 

An in situ flux experiment was undertaken in April, 1986. Lucite chambers 
were placed on the sediment surface and then pushed downward a few cm into the 
sediment. A battery operated stirring system was used to keep chamber waters mixed. 
Overlying water samples within the chamber were collected for SC0 2 concentration 
and isotope analyses shortly after emplacement and again at the end of the experiment 
prior to the next low tide ebullitive event. Cores were taken adjacent to the chambers 
and analyzed for EC0 2 concentration and isotope profiles. 

Flux experiments were performed in the laboratory in February 1987 and 
March 1988. Lucite box core chambers were used to collect an undisturbed section of 
sediment and overlying water. The chambers were sealed from below and transported 
in seawater back to the laboratory and placed in a seawater circulation tank to 
maintain near ambient temperatures. Overlying water samples were taken immediately 
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after the chambers were stabilized and sampled again twice during the following day. 
Chambers were circulated at a rate of approximately 15 ml/minute by circulating the 
overlying water with a peristaltic pump. At the end of the experiment the chambers 
were subcored and the porewater was extracted for SC0 2 concentration and isotopic 
analyses. During the March, 1988 experiment, overlying water and sediment samples 
were also collected for CH^ concentrations and isotopic analyses. 

RESULTS 

Particulate Organic Carbon 

POC content decreases exponentially with values ranging from 4% at the 
surface to less than 3% at depth (Fig. 2.1a). The 6 13 C values of the POC from this 
study as well as profiles from Blair et al. (1987) and Blair and Carter (1992) (Fig 
2.1b) behave conservatively as a function of depth, indicating that the fraction 
remineralized must have a 6* 3 C value similar to the site average of -19.08 ± 0.26 
(Blair and Carter, 1992). More formally, a mass balance calculation can be used to 
estimate the 6 13 C signal of the remineralized organic carbon: 

x/ , ( (%Sremi n ) M 10° -%C remin ) (Juried) >_ (2) 

og 0 ‘ Too 

where %G remin is the percent of the POC remineralized. Martens and Klump (1984) 
determined the percent of organic matter remineralized to be 28 ± 7%, based on five 
POC profiles. The isotopic composition of the POC at the surface (8 Go) is 
-19.1 ± 0.2 and at depth (6G buried ) is -19.2 ± 0.1 (Fig 2.1b). Solving Eqn. (2) for 
jG rem in gives an isotopic signature of the remineralized fraction of the organic matter 
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of -19.2 ± 0.2 per mil. 

SC0 2 Diffusive Flux 

The results of the SC0 2 flux chamber experiments (Table 2.1) were integrated 
with previous flux measurements (Martens and Klump, 1984) to yield an estimated 
annual flux of 30.3 ± 4.8 moles-m' 2 -yr _1 . The error estimate of ± 16% is based on 
the interannual variability of March flux measurements in 1977 , 1978 (Martens and 
Klump, 1984), and 1988 (this study). The SC0 2 diffusive flux across the sediment- 
water interface represents 80% of the remineralized organic carbon that is recycled to 
the overlying water. This value is consistent with the 84 ±18% measured by 
Martens and Klump (1984). The isotopic signal of the annual diffusive EC0 2 flux is 
estimated to be -15.5 ± 1.2 per mil based on the results of directly measured fluxes 
(Table 2.1). 

Fluxes and isotopic signatures were calculated from pore water gradients for 
the three flux experiments to compare to the measured flux values (Table 2.1). 
Calculated fluxes were determined using Fick's 1st Law of diffusion modified for 
chemical transport in sediments, 

J diffusive' "♦o ^ ^ 0 D s ^ 3 ^ 

where: 3 3 

<J) Q = sediment porosity at sediment surface (0.926 cm pw /cm* wet sed.) 
Martens and Klump, 1980) 

D s = sediment diffusion coefficient (cm 2 / sec) (D 0 taken from Li and Gregory 
s (1974) and corrected for sediment tortuosity (Berner 1980)) 

(3C/3z) 0 = linear gradient of SC0 2 , E 12 C0 2 or E 13 C0 2 profiles at the sediment 
water interface (mM/cm) 
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Calculated and measured fluxes and their isotopic signatures are compared in Table 
2.1. The calculated isotopic signal of the March EC0 2 flux was also determined by 
solving for the individual E 12 C0 2 and E 13 C0 2 concentrations and then calculating the 
distribution and isotopic signature of each of the following species based on sediment 
temperature, salinity, pressure and pH: 12 C0 2 , 13 C0 2 , H 12 C0 3 , H C0 3 , C0 3 , 

13 C0 3 = ; Deines et al., 1974; Stumm and Morgan, 1981) (Table 2.2; See appendix 
for description). pH data were taken from Chanton (1985). Individual gradients for 
each of these species were calculated and the corresponding diffusion coefficients were 
used in Eqn. 3 (Unver and Himmelblau, 1964; Ii and Gregory, 1974; Friedman and 
O'Neil, 1977; O’Leary, 1984). The resulting isotopic signal of EC0 2 flux was 
compared to a flux based only on the gradients of E l2 C0 2 and E l3 C0 2 (Table 2.2). 
The individual species were calculated because previous studies had found significant 
diffusion of C0 3 = into the sediment despite the net flux of SC0 2 out of the sediment 
(Sayles and Curry, 1988; McNichol et al., 1991). There was no difference in the 
calculated isotope signal of the EC0 2 flux using either approach and so individual 
species gradient calculations were not applied to the February and April fluxes. 
Application of a range of pH values (6-8) for the 0-1 cm interval shifted the resulting 

isotopic signal by less than 1 per mil. 

The isotope signature of the calculated EC0 2 flux was consistently 6 to 7 per 
mil enriched in 13 C compared to the measured values (Table 2.1). The measured 
fluxes were slightly larger than the calculated fluxes except in the April in silU 
experiment in which the flux was calculated using porewater gradients from a core 
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near the chamber. 


Burial of SC0 2 

Burial of SC0 2 accounts for a flux of 6.6 ± 0.5 moles-m* 2 -yr 1 with an 
isotopic signature of +8.3 ±1.5 per mil. The burial flux was determined using an 
average of EC0 2 concentrations measured at depths between 35 and 40 cm depth (Fig. 
2.2) and applying the following equation: 

(4) 

where 

= sediment porosity at depth (0.85 cm^ ^Jcvc? wet ^ ; Chanton, 1985) 

(ji m = sediment accumulation rate at depth (10 cm/yr: Chanton et al., 1983), 
C„ = average concentration of SC0 2 at 35 to 40cm (mM/cm). 

The burial depth of 40 cm was chosen because the SC0 2 concentration profiles are 

nearly asymptotic and greater than 95% of the remineralization occurs above this 

horizon (Martens and Klump, 1984). The isotopic signal of the buried EC0 2 of +8.3 

+ 1.5 per mil was determined by averaging the isotopic values measured at depths of 

35 to 40 cm. 

CH 4 Bubble Flux 

Methane ebullitive fluxes were measured previously (Martens and Chanton, 
1989; Martens et al., 1986). The annual ebullitive flux is calculated to be 6.4 ± 0.8 
moles-m^-yr' 1 based on fluxes measured from 1976 to 1986 (Table 2.3). The 
isotopic signature of the ebullitive flux of -60.0 ± 1.2 is based on flux-weighted 
average monthly methane measurements (Martens et al., 1986; Fig. 2.3). The isotope 
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measurements (Fig. 2.3) show strong seasonality, with relatively 13 C-ennched 
methane consistently released during the period of peak production (July- August). 

C0 2 Bubble Flux 

Ebullition strips a portion of the dissolved C0 2 from the porewater. The C0 2 
ebullitive flux, 0.13 ± 0.05 moles-m^-yr' 1 has an annual 6 13 C value of -8.5 ± 1.4 
per mil (Martens et al., 1986). The more negative C0 2 6 l3 C values were measured 
in the summer months corresponding to the period of highest remineralization rates 
and the most depleted SC0 2 values measured in the sediment (Fig. 2.2; Martens et 

al., 1986). 

CH 4 Diffusive Flux 

The annual methane diffusive flux of 0.85 ±0.8 moles-m 2 -yr 1 was taken 
from Martens and Klump (1980) and is based on pore water methane profiles and 
saturation calculations. The measured isotope signature of the methane flux is 
estimated to be -50.5 ± 0.3 per mil, based on the two flux chamber measurements 
from this study. The diffusive signal is 9.5 per mil enriched in 13 C relative to the 
annual methane ebullient isotope signal (Martens et al., 1986; this study) and is 4 per 
mil enriched compared to the diffusive flux calculated from pore water gradients of 

methane for March (-54.0 ± 0.3 per mil; Table 2.1). 

No attempt to determine the seasonality of this signature was made, however, a 
sensitivity of this flux to the overall isotopic signature of the remineralized fraction 
was performed. Assuming that either all or none of the diffusive flux was oxidized 
(the probable cause for the isotopic signature that is different from the methane 


produced in the sediment) resulted in less than a 0.5 per mil shift in the isotopic 
signature of the remineralized carbon. 

CH 4 Burial Flux 

The burial rate of methane has been calculated to be 0.14 ± 0.02 moles-m - 
yr' 1 (Martens et al., 1993). The 6 13 C value of the buried CH 4 (-59.9 ± 1.9) was 
determined by averaging the isotopic signal of the sedimentary methane at 40 cm (Fig. 
2.4). This signal is not significantly different than the methane bubble flux value of 
-60.0 ± 1.2 per mil, indicating that the ebullitive and burial processes are not 
influenced by oxidation or transport isotope effects (Martens et al., 1986; Chanton and 
Dacey, 1991). 

Solid Phase Inorganic Carbon 

Ca + + concentration profiles increase from overlying water values in the upper 
few centimeters and then decrease with depth (Fig. 2.5a). The rates of dissolution and 
precipitation of particulate inorganic carbon (PIC) were estimated from these 
porewater Ca + + profiles by assuming that an increase in the Ca + + concentration 
profile was caused by dissolution and a decrease in concentration was caused by 
precipitation of PIC. The Ca + + porewater linear concentration gradients were used 
in the 0-5 and 5-25 cm intervals to estimate the dissolution and precipitation rate of 
CaC0 3 respectively. The bulk sediment diffusivity for Ca + + (Dg) was determined 
by correcting the diffusion coefficients for Ca + + measured in seawater (Li and 
Gregory; 1974) for the given temperature and sediment tortuosity (Berner, 1980). 

The isotopic composition of the PIC that was assumed to have dissolved in the 
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0-5cm depth interval was determined by averaging the measured isotopic composition 
of the PIC within the dissolution zone (Fig. 2.5b). The isotopic composition of the 
assumed precipitated carbonate of -3.8 ± 3.8 per mil is based on the average isotopic 
value of the DIC pool in the 5-25 cm depth interval from the same months as the 
measured Ca + + profiles (see Fig. 2.2). The use of the Ca + + profiles to determine 
the PIC flux implies that CaC0 3 is the dominant species involved in the dissolution 
and precipitation processes. To a first approximation, this appears to be true since 
Mg ++ concentration profiles showed no significant trends and Sr ++ concentrations, 
which mimicked the Ca + + profiles, were an order of magnitude lower in 
concentration than the Ca ++ values (see Appendices for Mg ++ and Sr ++ results). 
The decrease in Ca + + concentration at depth may be a result of precipitation after 
core recovery because of the oversaturated nature of the porewaters, thus the calculat- 
ed precipitation rate should be regarded as an upper limit. 

The Carbon Isotope Mass Balance 

The carbon isotope mass balance determined as described above is summarized 
in Table 2.4 and shown schematically in Fig. 2.6. The total flux (J tota i) is defined as 
the sum of the buried, bubble and diffusive fluxes of SC0 2 and CH 4 . To determine 
the isotopic composition of the portion of J tota i that results directly from 
theremineralization of the POC (J remin ), the contributions of carbonate dissolution and 

precipitation need to be accounted for: 


J to tal~ Jr emit ! + ^diaa + ^procip 


( 5 ) 


{ ^ total ) t ® total^ = (Jiemin) + (^diss^ ^ dlse ) + ^piecip^ ^ precip ) (®) 


Using eqns. 5 and 6, the isotopic composition of the remineralized carbon is estimated 
to be -18.9 ± 2.7 per mil, which is in excellent agreement with the remineral i zed 
organic matter signal determined from the POC profiles (-19.2 ± 0.2 per mil). 

The error associated with the estimated signature of the remineralized carbon 
pool of 2.7 per mil is primarily due to interannual variability in fluxes and the errors 
associated with the isotopic composition of the diffusive ECO2 flux. The annual SCO2 
flux makes up approximately 70% of the remineralized carbon fraction, and therefore, 
the budget is very sensitive to the isotopic value assigned to this flux. The internal 
consistency of the mass balance suggests that despite the small data sets used to 
determine the methane and ECO2 diffusive fluxes and the interannual variability, the 
values are realistic estimates of the isotopic signature of these fluxes. 

The isotopic composition of the diffusive CH4 flux was calculated from a 
single experiment conducted in March. There is undoubtedly some seasonality to this 
signal that is not included in the estimate used here, but as noted, the small size of this 
flux precludes it from strongly affecting the overall isotopic signature of the 
remineralized carbon. 

The isotopic composition of the dissolved organic carbon (DOC) flux has not 
been measured and consequently is not included in this budget. The diffusive and 
burial DOC fluxes should not significantly influence the overall isotope mass balance 
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because they represent only 5% of the remineralized carbon flux (Martens et al., 

1993), and the isotopic composition of DOC fractions in similar environments are 
typically within 2-3 per mil of the particulate organic carbon pool (Williams and 
Gordon, 1970; Nissenbaum et al., 1972; Brown et al., 1972; Orem et al., 1986; 

Alperin, 1988). 

DISCUSSION 

The successful isotopic mass balance of the isotopic signature of the sources 
and sinks of carbon verifies the carbon budget determined for Cape Lookout as well as 
the isotopic signatures of the different carbon fluxes at a site that experiences 
significant seasonal variations. More importandy, the isotope mass balance allows us 
to constrain further our interpretations of certain processes occurring in this system to 
a degree previously not possible. These processes are discussed below in the context 
of the balanced carbon isotope budget. 

Surficial processes influencing EC0 2 and CH 4 fluxes 

The diffusive fluxes of both EC0 2 and CH 4 and their associated isotopic 
signatures were direcdy measured in chamber experiments. Pore water analyses of the 
flux chamber sediments were made so that the fluxes and their isotopic signature could 
be estimated from the concentration gradients at the sediment-water interface using 
Fick's 1st Law (Eqn. 3). Our original intent was to use the monthly SC0 2 profiles 
(Fig. 2.2) to determine the annual 6 13 C value of the diffusive flux. The comparison 
of measured and calculated SC0 2 fluxes and their associated isotopic signatures 
revealed that the calculated isotopic compositions of the fluxes were consistently 6-7 
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per mil enriched in 13 C relative to the measured values (Table 2.1). The calculated 
diffusive flux of CH 4 was greater and depleted in 13 C relative to the measured CH 4 
flux (Table 2.1). Several possible mechanisms that could account for the 
discrepancies between the calculated and measured isotopic signatures of these fluxes 
are discussed below. 

In theory, differential diffusion of the inorganic carbon species, C0 2 , HC0 3 , 
and C0 3 = , because of pH gradients across the sediment-water interface, could 
influence the isotopic composition of the flux of SC0 2 and the sedimentary SC0 2 
pool. The apparent production of EC0 2 that is significantly enriched in 13 C relative 
to the organic carbon has been noted in coastal and deep-sea sediments and has been 
interpreted to result from the diffusion of 13 C-depleted C0 2 out of and relatively 13 C- 
enriched C0 3 = into the sediment (McNichol, 1986; Sayles and Curry, 1988; 

McNichol et al., 1991). Calculations of the individual fluxes and isotopic signals of 
C0 2 , HC0 3 ' and C0 3 = from the March flux experiment demonstrate that the 
magnitude of the C0 3 = flux is insufficient to alter significantly the overall isotopic 
EC0 2 signal and thus this mechanism can be ruled out for Cape Lookout sediments 
(Table 2.2). 

Differences in the diffusion coefficient of isotopically substituted species 
would, if they exist, contribute to apparent diffusion isotope effects. For example, the 
ratio of the diffusion coefficient for 12 C0 2 and 13 C0 2 in water ( 12 D/ l3 D) is 1.0007 ± 
0.0002 at 25 °C (O’Leary, 1984). This fractionation has been included in the 
calculations of the SC0 2 isotopic signal (Table 2.1). To our knowledge, the potential 
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isotope effects associated with the diffusion of HC0 3 \ C0 3 or CH 4 have not been 
measured. 

We consider methane oxidation in surficial sediments to be a source of C- 
depleted material to the diffusive EC0 2 flux. Methane oxidation could also explain 
the discrepancies in the measured and calculated CH 4 fluxes. Evidence for CH 4 
oxidation includes 14 CH 4 -tracer studies (M. Alperin, pers. comm.), and the C- 
depleted CH 4 in the winter surface sediments at Station A-l (Fig. 2.4). Methane 
oxidizing bacteria preferentially utilize 12 CH 4 (Silverman and Oyama, 1968; Barker 
and Fritz, 1981; Coleman et al., 1981; Whiticar and Faber, 1986; King et al., 1989). 
The fractionation factor ( 12 k/ 13 k) for aerobic oxidation of CH 4 ranges from 1.005 to 
1.031 in culture (See Table 1 in Whiticar and Faber, 1986). An estimate of the 
fractionation factor for methane oxidation in the surficial sediments of CLB can be 
determined if it is assumed that the difference between the measured and calculated 
methane isotope signal of the March flux experiment is entirely the result of methane 

oxidation: 

&caic = ^moaa^~ F ox) ) + ^ 

where F ox is the fraction of the potential diffusive flux of methane that has been 
oxidized, and 6 calc and 6 meas are the isotopic compositions of the calculated and 
measured fluxes (-54.0 ± 0.3 and -50.5 ± 0.3, respectively) (Table 2.1). F ox was 
estimated to be 0.32 using the equation: 

Accordingly, S ox is -61.5 per mil and the fractionation factor (a) can be estimated 
using the following relationship: 
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( 8 ) 



a = 


(5,™* + l0 3 ) 


(5^+lQ 3 ) 


(9) 


The calculated a is 1.012 ± 0.005, which falls within the range of both the culture 
studies mentioned above, and the range determined by Whiticar and Faber (1986) for 
methane oxidation in natural systems (a = 1.002 to 1.014). 

If it is assumed that the difference between measured and calculated CH 4 fluxes 
is due to CH 4 oxidation alone, then only 14 p. moles- m ^-yr 1 of CO 2 are added to the 
SCO 2 flux with an isotopic signal of -61.5 per mil (Table 2.1). This small addition of 
13 C-depleted C0 2 shifts the calculated SC0 2 flux by less than one per mil, and thus 
does not fully explain the discrepancy between the measured and calculated isotopic 
values of the EC0 2 flux. 

Diagenetic Isotope Effects 

The carbon isotope mass balance allows us to address the effects of early 
diagenesis on the isotopic composition of the organic matter that is deposited at Cape 
Lookout Bight. For the purposes of this discussion, a diagenetic isotope effect is 
defined as any process that alters the isotopic composition of the organic carbon or 
creates isotopically distinct pools of carbon after the original organic matter is 
deposited in a sediment. Selective degradation of isotopically heterogeneous fractions, 
or kinetic isotope effects associated with various diagenetic processes could result in a 
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diagenetic isotope effect. The well known isotopic discrimination associated with 

biogenic methane formation is included in this definition. 

In principle. Cape Lookout Bight is an ideal location to investigate potential 
isotopic effects because of its nearly steady state annual input of organic matter. In 
addition, Cape Lookout Bight receives organic matter from a variety of isotopically 
distinct sources such as phytoplankton and seagrasses, as indicated by isotopic 
measurements, lignin and lipid analyses (Haddad, 1989; Blair and Carter, 1992). 
Isotopic differences apparently exist between compound classes at this site as well, as 
indicated by the average 6 13 C values of fatty acids (-22.1 ± 0.5), neutral lipids (-22.9 
± 0.3) and POC (-19.08 ± 0.26) fractions (Blair and Carter, 1992). Thus, given that 
nearly 30% of the organic matter delivered to site A-l is remineralized, a diagenetic 
isotope effect might be expected, if any of these organic fractions were selectively 
degraded. Surprisingly, the downcore profiles of the POC 6 13 C values and the 
isotope mass balance of the SC0 2 and CH 4 fluxes provide no evidence for selective 
degradation or preservation of isotopically distinct organic pools. McArthur (1989), 
Cronin and Morris (1982), Reimers and Suess (1983) and others, also found no 
indication of an isotopic diagenetic effect in the organic carbon record from a variety 
of sediment types and ages. This is in contrast to observations made in marine 
sediments of Skan Bay (Alperin, 1988) where changes in the isotopic composition of 
the organic matter downcore could be attributed to different rates of remineralization 
of the dominant organic matter sources, and Mangrove Lake, Bermuda where a C 
enrichment downcore in the sediments was attributed to selective preservation (Spiker 
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and Hatcher, 1984). 

The absence of any expressed isotope effect associated with the input and early 
diagenesis of organic matter at CLB is surprising since the inputs are isotopically 
distinct. The analyses of both the deposited and buried fractions of the organic matter 
of Haddad (1989) and Martens et al. (1993) indicated that the organic matter that was 
remineralized was most similar to an algal/bacterial source based on the C/N ratios 
and the loss of biochemical components downcore. The carbon isotope mass balance 
can be explained most simply by a remineralized fraction that is a mixture of two 
sources, predominantly algal/bacterial and a vascular plant component. 

The large isotopic difference observed between the POC and CH 4 indicates that 
isotope fractionation occurs during some step or steps in the production of methane 
from the breakdown of organic matter. Large isotope effects are associated with the 
reduction of C0 2 and the dissimilation of acetate to form CH 4 at this site (Blair and 
Carter, 1992; Blair et al., 1993). The carbon isotope budget should allow us to 
determine if C0 2 reduction and acetate dissimilation are masking other diagenetic 
isotope effects occurring during the initial breakdown of organic matter at Cape 
Lookout. 

The carbon isotope mass balance, when considered in light of the downcore 
6 l3 C acetate measurements of Blair and Carter (1992) indicate that there is little or no 
fractionation in the breakdown of high molecular weight compounds (biopolymers) to 
form the low molecular weight biomonomers. If there was isotopic fractionation in 
the initial breakdown of the organic matter, one would expect to see 6 C acetate 
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measurements in the sulfate zone that were very different from the 6 C POC 
measurements downcore, and the isotopic composition of the flux weighted sum of 
organic carbon remineralization (6 total ) would be isotopically different than the organic 
matter. Isotope discrimination might occur during the fermentation of lower weight 
organic matter to form acetate and other intermediates, however its signal may not be 
expressed if most of the organic matter flows through the acetate pathway (Blair and 
Carter, 1992). An acetate 6 13 C measurement from the sulfate reduction zone of Cape 
Lookout Bight sediments of -17.6 per mil indicates that the expressed isotope effect is 
probably less than 2 per mil during the oxidation of acetate (Blair and Carter, 1992). 
These observations suggest that the isotopic fractionation of the carbon during the 
formation of methane via C0 2 reduction and acetate dissimilation occurs in the final 
steps of CH 4 production and little or no fractionation is expressed in the breakdown of 
higher molecular weight compounds to form low molecular weight compounds such as 
acetate. This conclusion is in contrast to the less well-constrained speculations about 
sedimentary processes that would cause intramolecular isotope effects such as bond 
rupture (Macko, 1992), polymerization, and decarboxylation (Galimov, 1980) and thus 
alter the isotopic composition of the residual organic matter. 


SUMMARY AND CONCLUSIONS 


The carbon isotopic signal of the remineralized organic matter in Cape Lookout 
Bight sediments is -18.9 ± 2.7 per mil based on a mass balance of the remineralized 
carbon fractions. This estimate is indistinguishable from the calculated 5 l3 C signal of 
the remineralized organic matter determined from vertically uniform POC profiles of 
-19.2 ± 0.2 per mil. The agreement confirms that the major sources and sinks of 
carbon and their isotopic signals have been quantified. The comparison of directly 
measured and calculated CH 4 and SCO 2 isotope fluxes suggest that a surficial process 
is altering the isotopic composition of these species diffusing out of the sediment. 
Methane oxidation has altered the isotopic compositions of both the E CO 2 and CH 4 
leaving the sediment, but can account for little of the discrepancy in calculated versus 
measured 6 13 C values of the EC0 2 fluxes. 

The isotope mass balance revealed little variation in the isotopic composition of 
the organic matter being deposited, remineralized, or buried. This is in contrast to 
numerous recent studies that have focused on organic matter source variations, 
individual biochemical components and their selective loss or preservation, and also 
intramolecular isotope effects to explain downcore changes in 6 l3 C values of POC 
(Galimov, 1980; Spiker and Hatcher, 1984; 1987; Benner et al., 1987; Alberts et al., 
1988; Hayes et al.,1990; McArthur et al., 1992). The absence of an expressed 
diagenetic isotope effect in a system with a mixture of organic matter sources suggests 
that the processes that have been speculated to cause the isotopic variations in 
sedimentary organic matter profiles may not typically be expressed in the sedimentary 
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record, and the variations seen may more often be source related (Dean et al., 1986). 

A carbon isotope mass balance was achieved by seasonal measurements of the 
dominant components of the carbon cycle in this system, highlighting the importance 
of considering seasonality when determining an isotopic signal for a specific 
environment. This is particularly important in coastal systems where the rates of 
remineralization processes and the isotopic signals of the various fluxes should vary 
seasonally. 
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Fig. 2.1 Cape Lookout Bight particulate organic carbon concentration (A) and 
6 13 C values as a function of depth in the seabed. (B) POC profiles are 
from Chanton et al. (1983); Blair et al. (1987), Blair and Carter (1992) 
and this study. 
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Fig. 2.2 SC0 2 concentration and 6 13 C profiles from Cape Lookout Bight. 

Samples were collected in 1986 except where noted. Individual 
symbols represent separate cores collected on the same sampling date. 
The open and filled circles for the May 1986 profiles represent 
duplicate samples from the same core. Note differences in 
concentration, depth, and 6 13 C scales for March profiles. 
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Fig. 2.3 Methane flux and & 13 C bubble data from Cape Lookout. The 

ebullitive fluxes from mid October to May are essentially zero. Flux 
and 1983-1984 6 13 C bubble data are from Martens et al. (1986). 6 C 
bubble data for 1986 are from this study. 
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Fig. 2.4 Sedimentary CH 4 S 13 C profiles from Cape Lookout Bight. 
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Fig. 2.5 Cape Lookout Bight porewater calcium profiles from 1990-91 (A) and 
partic ulate inorganic carbon concentration (PIC) and 5 C profiles (B). 
The dotted line at 5 cm indicates the maximum depth assumed for net 
dissolution and 5 to 25 cm is considered the zone of precipitation. 
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Fig. 2.6 The carbon flow and isotopic composition of the major fluxes at CLB. 

The double arrows for methane and COj represent diffusion and bubble 
ebullition. Carbon fluxes are given in parentheses. 
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diffusion coefficient taken from Alperin (1988) 
2 d for temperature and porosity. 




Table 2.2 Concentration and <S 13 C of ZC0 2 and individual dissolved inorganic species 
in flux experiments for March 1988. 
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Table 2.3 


CH 4 Bubble Flux from Martens et al. (1986) and Martens and 
Chanton (1989). 


(moles-m^-month' 1 ) 


Month 

1976 

1977 

1978 

1985 

1986 

Average 

s.d. 

May 




0.08 

0.54 

0.31 

0.32 

June 


0.98 


1.23 


1.10 

0.18 

July 


0.31 


2.21 


1.26 

1.35 

Aug 

0.70 

1.71 

1.57 

2.94 

2.02 

1.79 

0.81 

Sept 

0.46 


0.86 

2.75 


1.36 

1.22 

Oct 

0.71 


0.16 

0.89 


0.59 

0.38 



Sum (mol-m' 2 - 

•yr' 1 ): 


6.4 

2.1* 


*Error was determined as: 


o=-/ETs7d7P 
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Table 2.4 Carbon isotope mass balance for Cape Lookout, Station A-l. 


moles- m' 2 -yr _1 fi l3 C 


ch 4 



Buried 

0.14(0.02) 

-59.9(1.9) 

Bubble 

6.4(2. 1) 

-60.0(1.2) 

Diffusive 

0.85(0.8) 

-50.5(0.3) 

sco 2 



Buried 

6.5(0.5) 

+8.3(1.5) 

Bubble 

0.13(0.05) 

-8.5(1. 4) 

Diffusive 

30.3(4.8) 

-15.5(1.2) 

^total 

44.3(12.6) 

- 19 . 2 ( 2 . 1 )' 

PIC 



Dissolved 

0.27(0.4) 

- 0.14(0.2) 

Precipitated 

-1. 1(0.3) 

-3.8(3. 8) 

J • 2 

•'remm 

45.13(12.6) 

-18.9(2.7) 


lr The estimated fi 13 C of the flux (J total ) was determined using the following equation: 


5 13 C(J totaJ ) = 


2 1(5^) (J,)] 
2(^1 


where i = the diffusive, ebullitive and burial fluxes of CH 4 and 2C0 2 . 
Calculated using Eqn. 5. 
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